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ABSTRACT 

The electronic structure and optical properties of metallic nanoshells are investigated using a jellium model and 
the Time Dependent Local Density Approximation (TDLDA). An efficient numerical implementation enables 
applications to nanoshells of realistic size with up to a million electrons. We demonstrate how a frequency 
dependent background polarizability of the jellium shell can be included in the TDLDA formalism. The energies of 
the plasmon resonances are calculated for nanoshells of different sizes and with different dielectric cores, dielectric 
embedding media, and dielectric shell backgrounds. The plasmon energies are found to be in good agreement 
with the results from classical Mie scattering theory using a Drude dielectric function. A comparison with 
experimental data shows excellent agreement between theory and the measured frequency dependent absorption 
spectra. 
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1. INTRODUCTION 

The optical properties of metallic nanoparticles are of considerable fundamental and technological interest.^ A 
particularly interesting nanoparticle is the metallic nanoshell,^ which consists of a metallic layer grown over a 
solid dielectric core. The plasmon frequency, which determines the nanoshell optical properties, can be tuned 
over a wide spectral range by simply varying the ratio of inner to outer diameter of the shell. The tunability of the 
optical properties of metallic nanoshells enables several important applications such as resonant photooxidation 
inhibitors,^ optical triggers for drug delivery implants,^' ^ environmental sensors,^ and Raman sensors.^ 

Although classical Mie scattering seems to account well for the optical properties of metallic nanoshells, there 
is clearly need for a more microscopic description. A microscopic understanding of the electronic properties of 
nanoshells is a prerequisite for a systematic and rational modification of their electronic and optical properties. 
We have recently started the development of a fully quantum mechanical model for the electronic properties of 
metallic nanoshell using density functional theory. ^^^'^ In this paper we will review the general features of our 
approach and present some new results. The emphasis will be on the effects of dielectric backgrounds on the 
electronic structure and optical properties of metallic nanoshells. 

2. ELECTRONIC STRUCTURE 

TDLDA is a well established method,^* which allows surprisingly accurate calculations of the electronic and 
optical properties of nanoparticles. The method was extensively applied to small metal clusters and proven to 
be in excellent agreement with more exact calculations and experiment. The method was recently imple- 
mented to shell geometries and proven to be in good agreement with the existent classical and semi-classical 
approaches. ^""^'^ TDLDA involves two steps. The first step consists of self-consistent calculation of the electronic 
structure and the independent electron response functions of the nanoparticle. The second step consists of solv- 
ing self-consistently the RPA equation for the screened response function of the system. Our implementation of 
the first step to shell geometries has been discussed in details previously, ^"^^ and therefore only the main points 
will be presented here. 
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Figure 1. The effect of different dielectrics on the electronic structure of a (70,100) a.u. gold nanoshell. The curves in 
each panel represent, from the bottom and np, the effective potential, the Coulomb potential (the values are on the left 
axis) and the electron density (the values on the right cixis). Different panels correspond to: a) £c = 1, = 1 and sj = 1; 
b) EC = 5, eb = 2 and ej = 1; c) ec = 1, eb = 1 and ej = 8; d) ec = 5, eb = 2 and sj = 8. For comparison, in panels 
b)-d) we have included with dashed lines the electronic structure when no dielectrics are present, i.e. ec = £e = £j = 1. 



The electronic structure of the nictallic nanoshells is determined from the Kohn-Sham equations 

(^-^^ + Vext+VH +Vxc[n\^ ^i=£i(t>i, (1) 

where 

n(f)=^(l + e'5(^'-''))"'|</..(a-)^ (2) 

i 

The chemical potential is fixed by the condition: 

7V = ^(l + e^(^*-''))"\ (3) 

i 

N being the number of conduction electrons. The term Vext represents a uniform background potential inside 
the shell. Its value is fixed such that the converged calculations lead to a ionization potential equal to 5.4 eV, 
appropriate for gold. To avoid problems associated with the large degeneracies in spherically symmetric systems, 
we perform the calculations at low but finite temperature. The self-consistency is achieved by employing the usual 
iterative method. In the jellium approximation, the system becomes spherically symmetric and consequently the 
Kohn-Sham orbitals must be of the form r~^uik (r) Yim {f). In this case the Kohn-Sham equations are reduced 



to a set of coupled radial Schrocdingcr equations: 

/ 1 d2 la + l) , A , , 

\2d^ + ^'^^^t^J J "'^^'■^ " ^ikuikir), (4) 

where 2 1 21 1 

Even with the spherical symmetry, in order to make contact with the experiment, one still needs to calculate 
a very large number of orbitals. For example, in our largest simulation performed so far, the angular quantum 
number I was as high as 300 and the radial quantum number k of the discrete states was as high as 18. The 
dimensions of this nanoshell were (16,19) nm and, from an experimental point of view, this will be considered 
a very small nanoshell. The success of such large calculations clearly depends on how fast each iteration is 
performed. 




Figure 2. The electronic structure of three different nanoshells which were experimentally fabricated and characterized 
by optical measurements. In each panel, the curves represent, from bottom and up, the effective and Coulomb potential 
(left axis) and the electron density (right axis). The dimensions of the nanoshells are: (4.1,5.1) nm (panel a), (8.6,9.9) 
nm (panel b) and (13.1,14.9) nm (panel c). 



The radial Kohn-Sham orbitals and energies are calculated using the shooting method. This method consists 
of finding two, generally independent, solutions of the radial Kohn-Sham equation corresponding to an arbitrary 
energy. One solution is regular in the origin and the other is regular at infinity.^ These solutions are found by 
direct integration of the radial Kohn-Sham equation using the fourth order Rungc-Kutta method. The eigenvalues 
correspond to those particular values of the energy where the two solution become linearly dependent, i.e. their 
Wronskian is zero. Wc notice that this approach provides the eigenvalues and the orbitals in the same time. To 
find all eigenvalues, one will have to sweep the energy from the bottom of the effective potential to zero and 
record the points where the Wronskian changes sign. If there were no relation between different eigenvalues this 
will represent a prohibitively time consuming step. Fortunately, for a shell geometry, the radial Schrodinger 
equations corresponding to consecutive angular quantum numbers I have very close eigenvalues and £;+i,/c is 
always above si^k- This means we can reduce the sweeping interval to small intervals around e/.fe. A: = 0, 1, 2, 
in order to find the eigenvalues corresponding to I + 1. Thus, all the I > orbitals can be calculated very 
efficiently. This is the reason why the shell geometry allows calculations with such large number of particles. 
Suppose we keep the thickness of the shell d constant and we increase the outer radius b of the shell. In this case 
the maximum value of the radial quantum number k remains constant while the number of particles increases 
approximately as = 3db^/rl. This means that, as the overall size of the shell is increased, only orbitals with 
higher angular quantum number appear and, as we argued above, they can be calculated very efficiently. 

Another difficult problem wc encountered in our large simulations is the convergence of the iterative process. 
It was already pointed out in the literature that the iterative methods in general have convergence problem for 
large many-body systems. ^'^ Our studies on the thermodynamic limit of the Kohn-Sham and Hartree equations 
revealed that this convergence problem is due to the long range of the Coulomb interaction. • ' • ^ The solution to 
this problem is to replace the Coulomb interaction with a screened Coulomb interaction in the Hartree potential: 

/g-|x-i7|/A 
|^^_ ^1 [n (y) - no\dy. (6) 

The iteration converges fast for a small screening lengths A and suppose the result is nx. The next step is to 
increase the screening length A by a finite amount and start again the iteration process with n\ as the initial 
input. This is repeated mitil the results stabilizes. In the present calculations, the screening length was increased 
by 2 a.u. and, for each A, 8 iterations were needed to achieve the convergence. The first four decimals of the 
eigenvalues become independent of A for A > 40 a.u.. 

Experimentally, the metallic shells are grown over a solid dielectric, which so far was either AU2S or silica, 
and they arc usually suspended in solution during measurements. Also, the gold ion cores arc highly polarizablc 
due to the completely filled d-band. Consequently, the conduction electrons are embedded in a highly polarizable 
medium. The effect of these dielectric media can be included in the electronic structure calculations as follows. 
Let us consider a nonuniform, possibly frequency dependent dielectric, 

{El = £c for r < ri 
£2 = sj for n <r <r2 (7) 
£3 = se for r2 < r, 

where ri and r2 represents the inner, respectively outer radius of the shell. This parameterization accounts 
for the presence of a dielectric core, polarizable jellium background (metal cores) and a dielectric embedding 
medium. The Hartree term, 
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which esscintially represents the cilcctrostatic potential generated by the conduction electrons and the positive ion 
cores, has to be modified such that it satisfies the usual boundary conditions at the interface between different 
dielectric media:^^ 

VHir)/ec + ^c for |r| < ri 
vnir) = { VH{r)/ej + $j for ri <\r\< r2 (9) 
VH{r^/sE for r2 <\r\, 



The constants are found from the requirement that the potential is continuous and they are given by; 

r $j = {I/ee - l/sj) Vh (r2) 
\ $c? = $j + (l/£j - l/ec) Vh (ri) 

In Fig. 1, we show the effect of the dielectric media on the electronic structure of the nanoshell. The values 
for the dielectric constants were chosen as ec = 5, = 2 and £j = 8 which are representative for realistic 
situations. The influence of ec and ee is relatively weak. The strongest influence on the electronic structure 
comes from e j which strongly reduces the dipole fields near the surfaces. This leads to a lowering of the effective 
potential in the nanoshell. The reduced workfunction results in an increased electron spill out at the surfaces 
of the shell. The effect of the d-electrons is slightly reduced when the core and embedding medium are added. 
Although the dielectric effects on the electronic structure are visible, these variations alone have little influence 
on the optical response of the nanoparticles.^° This is in contrast to what have been observed experimentally, 
i.e. large shifts in the plasmon frequencies due to the presence of different dielectric media. As will be argued 
in the next section, the explanation is that the system responds differently to an external electric field in the 
presence of the dielectrics which screen the excitation. 

Figure 2 illustrates the converged electronic structure of three nanoshells which were fabricated and charac- 
terized by optical measurements. The electronic structure shown in this figure will be later used to calculate 
the optical response of these nanoshells and make a direct comparison with the experiment. The dimensions 
of these nanoshells arc (4.1,5.1) nm, (8.6,9.9) nm and (13.1,14.9) nm. The largest nanoparticle contains about 
2.5 X 10^ conduction electrons. The nanoshells were suspended in an aqueous solution (e^ = 1.78) and their core 
consisted of AU2S (ec = 5.4). The value of the jellium dielectric constant was fixed to ej — 8, determined from 
a fit of the experimentally measured gold dielectric function with a Drude like expression. 

The independent electron response function is calculated directly from the electronic structure. We arc 
interested here only in the dipolar term {I = 1) of the response function, which is calculated using the following 
expansion: 

til=i{r,r';u) = ^ \t ^ X] (gfcii) ^fcii (r) Mfcu (/) 

|ii-J2|=l k 

X [Gi^ (r, r'; Skh +ui + i6) + Gi^ (/, r; Skh - u) - iS)] , 
where C?;, / = 0, 1, represent the radial Green functions: 

("^ ^ + + ""'^^ (r) - - ij) Gi (r, r'; E + iS) = -5 (r - r') , (12) 

and Uki are the self-consistently calculated Kohn-Sham orbitals. The Green functions are calculated by direct 
integration of their specific equation. Eq. (11) provides the response function multiplied by r^r'^ which is more 
convenient to use as it will follow from the next section. 
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3. OPTICAL RESPONSE 

We will focus in this section on deriving the RPA equations in the presence of the dielectric medium defined 
by Eq. (7). The physical picture can be described as follow. When the whole system is placed in an external 
field, the conduction electrons will try to screen the excitation field. The dielectrics will polarize when placed in 
the external field, screening at their turn the excitation field. The conduction electrons will screen not only the 

original excitation but also the electric fields produced by the polarized dielectrics. To quantify these processes, 
we start from the TDLDA expression of the screening charge induced by an external electric field £'oe~*'^* e^, 

6n{f, w)= [ rfr"'n(0)(f, f '; w) [5v,, {?', w) + <l>c (r ', w)] , (13) 



where n^"-* represents the response function of the independent electrons. The quantity 6vxc represents the 
variation of the exchange-correlation potential due to the screening charge and represents the Coulomb 
potential generated by the screening charges Sn and including the excitation field. In the quasistatic limit, which 
is generally accepted to be valid for particles smaller than 40 nm,^^' the electrtatic potential 0c must satisfy 
the Poisson equation with the appropriate boundary conditions: 



~ dii 



(14) 



(l)c (r, w) -fEo for |f| oo. 
The solution can be written in the following form: 



47r 



3 J r?.e(r') 



Sn 



(r, Lo) r'^dr' + J^i-i 2 ^^^^'^^i ) ^0 cos 9. 



(15) 



The surface charges o"^'^ arc found from the boundary condition as described in a previous application.^^ After 
the solution is plugged back into Eq. (13) and using the notation vi {r,r') = r^/r?,, V2 {r,r') = dv\/dr and 
a (r, Lo) = r'^Sn (r, a;), the RPA equation becomes: 



a (r, Lo) - I dr' n^°^ (r, r'; w) /"^^^ (n (r')) a (r', to) 
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dr' I dr"&^Hr,r';u;)^-^^^a{r",co) 



where 



EfHr,co)l'-^air',.)dr' 
dr'nl") (r,r';u;)(/.0;(r'), 



X(2)(r,a;) = 



2^e 



'^1 {r,^) + ^^^-^X2{r,u,). 



and 



Xi (r, w) = y dr' 11^°) (r, r'; a;) (r', n) . 
In these equations, (j)^ represents the unscreened response of the dielectrics to the electric field: 
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are the classical polarization charges and is just a dielectric factor. 
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Figure 3. The effect of different dielectrics on the optical absorption of the gold nanoshell shown in Fig. 1. Different 

panels refer to the following situations: a) ec ~ 5, Ee ~ ^ and ej — 1; b) ec = i, Ee = "2 and ej = 1; c) ec ~ 5, Ee ~ '2 
and ej = 1; d) ec ~ 1, £e = 1 and ej — 8; e) ec ~ 1, £e — 2 and ej = 8 (solid) and ec = 5, se = 1 and ej = 8 (dotted); 
f) EC = 5, eb = 2 and sj = 8. For comparison, in all panels we have included the absorbance of the gold nanoshell when 
no core or embedding medium is present and ej = 1 (dashed line). 



The integrals in the above equations are discretized and the integral equations are transformed in algebraic 
equations which arc solved by conventional methods.®' The photo absorption cross section of the nanoshell is 
related to its polarizability a(w) through 

i7a6s(w) = ^/e^-Im[a{w)], 
c 

where SI units were used. The frequency dependent polarizability is calculated from 

a{uj) = ^ j dr r^a (r, lo) . (23) 

where a(r, w) is the local polarizability derived from the RPA equation Eq. (17). The absorbance is defined by 

the ratio between the optical absorption cross section and Trrg. 

In Fig. 3, we show how the photoabsorption cross section of a (70,100) a.u. nanoshell is influenced by the 
presence of dielectric media. For each nanoshell, the photoabsorption spectrum is dominated by two peaks, 
corresponding to the symmetric and antisymmetric coupling between charges on the inner and outer surfaces of 
the shell. The nanoshell plasmons can be viewed as bonding and antibonding superpositions of a cavity plasmon 
located on the inner surface of the shell and a sphere surface plasmon on the outer surface of the shell. In the 

presence of dielectric backgrounds, the energy of the cavity plasmon is uc = ^b\/ ^ and the energy of 

the sphere plasmon is ujg = ^B\j ^^^^^ ^'^^ energy (long wavelength) symmetric nanoshell plasmon a;_ 
has a larger admixture of the sphere plasmon and the high energy (low wavelength) antisymmetric mode a;_|_ 
is primarily composed of the the cavity plasmon. The changes in the optical spectra in Fig. 3 caused by the 
different dielectric media can be simply understood from how the bare cavity and sphere plasmon shifts in the 
presence of dielectrics. The a;+ plasmon with its large cavity plasmon content is more strongly influenced by a 
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Figure 4. The optical response of a (60,90) a.u. gold nanoshell when the dielectric constant of the core and embedding 

medium arc varied over a range of 1 to 4. The dielectric constant of the jellium shell ej = 1. A logarithmic scale is used 
on the y=axis. Each panel, from the bottom curve and up, represents the optical response corresponding to a dielectric 
core of 1,2,3 and 4. Different panels corresponds to £b = 1 ( panel a), es = 2 (panel b), eb = 3 (panel c) and = 4 
(panel d). 

dielectric core sc than by the embedding medium. The ui- plasmon is more strongly influenced by the dielectric 
embedding medium ee than by the didc'ctric core. The presence of a dielectric jcUium background sj has a 
stronger influence on lv+ simply because the cavity plasmon depend more strongly on sj than the surface sphere 
plasmon. 

In Fig. 4 we calculate the optical absorption spectra of a (60,90) a.u. for a broad range of values of the 

core and embedding medium dielectric constant. The optical absorption is plotted in a logarithmic scale. The 
flgure show that in addition to the two collective resonances 0;+ and U- there is a background contribution to 
the optical absorption. This background absorption is caused by discrete single particle excitations. It can be 
seen that as the dielectric constant of the core and embedding medium is changed, only the collective plasmon 
resonances are influenced. The discrete single particle excitations remain essentially unchanged because the 
presence of dielectric media only weakly influences the electronic structure of the system. A calculation of the 
polarizability only including the changes in the electronic structure caused by the dielectrics has almost no effect 
on the optical properties of the nanoshell.^'' The large effect of a dielectric core or embedding medium on the 
plasmon resonances is due to the polarization of the dielectrics and their contribution to (pc in Eq. (13). 

The frequency of the plasmon resonances extracted from the optical absorption spectra shown in figure Fig. 4 
are in excellent agreement with the predictions of the classical Mie theory as it is illustrated in Fig. 5. In the 
classical calculations, the metallic phase is modeled by a pure Drude dielectric function corresponding to the 
same Tg = 3 as used in the TDLDA calculations. 

The fundamental property of the metallic nanoshells is their tunable plasmon resonances, which was demon- 
strated in a series of experiments where monodisperse nanoshells of different shell thickness and core sizes have 
been fabricated and characterized with optical measurements.^^ In the absence of any dielectrics, the classical 
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Figure 5. The TDLDA calculated plasmon energies (circles) compared with the prediction of the classical Mie theory 
(solid lines) for different core dielectrics and embedding media. 

Mie theory predicts the following expression for the plasmon energies: 

'^?± = ^[1±^V^1^4Z(ITI)^], (24) 

which depends only on the aspect ratio of the nanoshell x = ri/r2- The same expression can be derived using 
semi-classical approaches. We show in the following that our self-consistent calculations lead to the same aspect 
ratio dependence of the plasmon energies. For this, we fix the dielectric constants to 1 and wc apply the TDLDA 
formalism to nine nanoshells of the same thickness (17 a.u.) but with different overall diameters, ranging from 
85 to 680 a.u.. The aspect ratio of the nanoshell is varied in this way over a broad range of values starting from 
0.6 and up to 0.95. From the polarizability curves wc extract the wavelengths of the collective modes which 
are shown in Fig. 6 by open circles. The agreement between the TDLDA calculated plasmon wavelengths and 
Eq. (24) is almost perfect over the entire range of aspect ratios. The lower energy plasmon becomes strongly 
dependent on the aspect ratio for thin shells when the interaction between the cavity and sphere plasmon becomes 
very strong. This is the regime where the nanoshells become very useful in a variety of applications. 

4. COMPARISON WITH EXPERIMENTAL DATA 

In this section we apply the formalism presented in the last section and calculate the optical response of the 
three nanoshells presented in Fig. 2. The procedure for calculating an effective jellium background dielectric 
function consists of the following. The conduction electrons are assumed to contribute to the dielectric function 
of the metal through a Drude term: 

sd{oj) = l-LO%/{uj'^ + iju;). (25) 
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Figure 6. The TDLDA calculated plasmon wavelengths (circles) compared with the prediction of the classical Mie theory 
Eq. (24) (solid lines) for nanosliells of different aspect ratios. The TDLDA has been applied to a set of nine nanoshells 
of the same thickness (17 a.u.) and with overall diameters of 85, 97.2, 117.4, 136, 170, 226.6, 340, 453 and 680 a.u.. 

This assumption is supported by the excellent agreement between the TDLDA calculations and the classical 

Mic calculations when the metallic phase is modeled by a Drudc dielectric hmction. The width 7 should be the 
same as the artificial width that will be used in the TDLDA calculations.^® The experimentally measured bulk 
dielectric function of gold,^^ £^^^^(0;) can then be decomposed into:^^'^^ 

e^^nu:)^ej{u) + eD{oj)-l, (26) 

which provides a simple way of determining the jellium dielectric function £j{uj). The RPA equation Eq. (17) 
have been solved with the frequency dependent jellium dielectric function obtained by this procedure and the 
rcsiilts arc shown in Fig. 7 by solid lines. Wc also show a calculation using a frequency independent ej = 8 
obtained from a best fit of £j{uj) to a constant. The agreement with the experimentally measured optical 
absorption shown by circles is excellent for all three nanoshells. The experimental spectra show an additional 
peak for all three nanoshells which is located at the same wavelength of about 520 nm. This peak is generated 
by solid gold colloids also present in the solution during the experiment. The peak at larger wavelengths is the 
W- resonant mode and its position is different for the three nanoshells, which have different aspect ratios. The 
u!- mode in the TDLDA calculations is slightly red shifted relative to the experiment. 

In these experiments the nanoshells are not perfectly monodispersc? and this is due to slight variations in their 
geometry. The size distribution was evaluated to a Gaussian of a stairdard deviation of 11%, which can explain 
the difference in the widths of the plasmon resonances. To support this argument, we have calculated the optical 
response of nine nanoshells of sizes uniformly distributed in between ±11% of the standard size. The results 
are shown in Fig. 8. The optical absorption curves have been multiplied by the corresponding Gaussian weight. 
As one can see, a size distribution of 11% introduces a large broadening of the plasmons. This is mainly due 
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Figure 7. The optical response of three gold nanoshells with AU2S' core of geometries (4.1,5.1) nm (panel a), (8.6,9.9) nm 
(panel b) and (13.1,14.9) nm (panel c). The solid curves represent the TDLDA calculation when the full frequency 
dependent dielectric constant of gold is used. The dashed lines represent the TDLDA calculation when a Drude dielectric 
function with ej = 8 is used. The circles represent the experimentally measured optical response. The dotted lines 
represents the optical absorption of solid gold particles also present in the solution during the experiment. 

to the fact that the thickness of the nanoshells is very small compared with the overall size and, in this regime, 
the plasmon modes are very sensitive to small variations of the geometry as it was illustrated in Fig. 6. Also, 
the size distribution induces a blue shift because the height of the plasmon peaks increases with the size of the 
nanoshells. The figure clearly shows that the variation of the plasmon energies caused by the inhomogeneous 
size distribution can account for the experimental broadening of the plasmon resonance. 

5. CONCLUSIONS 

The optical properties of large nanoshells can be calculated using the TDLDA method. The results show that 
the plasmon energies are sensitive to geometry and can depend strongly on the presence of dielectric backgrounds 
such as a dielectric core, jellium background and embedding medium. The calculated energies of the plasmon 
resonances are in excellent agreement with classical Mie scattering. A comparison with experimental data show 
very good agreement between the theoretical and experimental absorption spectra. 
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Figure 8. The cflcct of a size distribution on the optical absorption of a solution of gold nanoshells. The figure illustrates 
the TDLDA optical absorption of different nanoshells of slightly different geometries, 5 = ±2.75% (dashed), 5 = ±5.5% 
(dotted), 5 = ±8.25% (dashed-dotted) and S = ±11% (dashed-dotted-dotted) relative to the geometry (4.1,5.1) nm 
(solid). Each curve has been multiplied by a Gaussian factor ae"'''^''' , where a = 11% as predicted from experiment 
and a was chosen equal to 0.5. The circles represent the experimentally measured optical absorption. 
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